Dynamical Mean Field Theory of the Antiferromagnetic Metal to Antiferromagnetic 

Insulator Transition 
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We study the antiferromagnetic metal to antiferromagnetic insulator using dynamical mean field 
theory and exact diagonalization methods. We find two qualitatively different behaviors depending 
on the degree of magnetic correlations. For strong correlations combined with magnetic frustration, 
the transition can be described in terms of a renormalized slater theory, with a continuous gap 
closure driven by the magnetism but strongly renormalized by correlations. For weak magnetic 
correlations, the transition is weakly first order. 



PACS numbers: 71.30.+h 

The correlation driven metal insulator transition 
(MIT) or Mott transition is one of the central problems 
of condensed matter physics. For generic band struc- 
tures, it occurs at a finite value of the interaction strength 
and constitutes a genuinely non perturbative phenom- 
ena. Recently, a great deal of progress has been made 
using the Dynamical Mean Field approach (DMFT), a 
method which becomes exact in the well defined limit 
of infinite lattice coordination However, all stud- 

ies so far, have been confined to the paramagnetic metal 
(PM) to paramagnetic insulator (PI) transition. In this 
letter, we study the transition from the antiferromag- 
netic metal (AM) to an antiferromagnetic insulator (AI) 
within DMFT. The motivation for this work is twofold. 
Experimentally, the interaction or pressure driven MIT 
in V 2 3 |[|] and NiS 2 - x Se x |]|§ take place between 
magnetically ordered states. While we are still quite far 
from a realistic modeling of these oxides, (in this work 
we only consider commensurate magnetic order and ig- 
nore orbital degeneracy and realistic band structure) we 
would like to take into account the effects of magnetic 
long range order present in these systems. On the the- 
oretical side, we would like to understand how magnetic 
correlations , which controls the scale at which the spin 
entropy is quenched, affect the MIT. 

We consider the Hamiltonian, 



(1) 



t\ and t 2 are nearest and next nearest neighbor hoppings 
respectively and U is the onsite coulomb repulsion. J is 
a ferromagnetic spin spin interaction that we add to the 
model in order to independently control the strength of 
the antiferromagnetic correlations of the system H . We 
parametrise the hoppings as if = (1 — a)t 2 and t%- = at 2 
where a is the degree of frustration. t 2 breaks the per- 
fect nesting characteristic of the bipartite lattice ^,^] and 
allows for a direct AM to AI transition. 

Within the dynamical mean field approximation, in the 
presence of magnetic order, the single particle Green's 



function on the A and B Neel sublattices has the follow- 
ing form 



G-' = 



(2) 



a = 1,2 refer to the spin up and down states, m the 
staggered magnetisation and fi a = (fi — ^ + (— l) a Jm). 
E are the self energies which are momentum indepen- 
dent within this approximation. The local Green's func- 
tions Gua- depend on the sublattices (A and B) and 
obey Gai = Gb2 = G%, Ga2 = Gbi = G 2 . Similarly, 
Eai = £_B2 = Si and XU2 = Sbi = £2- 

Within DMFT, these local self-energies and Greens 
functions are obtained from an Anderson impurity model 



(3) 



describing an impurity electron / hybridising with a bath 
of conduction electrons d via the spin-dependent hybridi- 
sation functions 



A a (iuj n ) = ^2 ' 



(4) 



where w n are the Matsubara frequencies. The A CT obey 
self consistency conditions which can be expressed in 
terms of G± , G 2 and the non- interacting density of states 
We restrict ourselves to the Bethe lattice in the 
following. The non-interacting density of states on the 
Bethe lattice is given by D(e) — -^jyzVD 2 — e 2 with 
D = V2t being the half bandwidth. The self-consistency 
conditions now have a simple form 



t 



ti 



A cr (iw„) = -^G- a (iuj n ) + -^G a {iuj n ) 



(5) 



We consider only the half-filled case, which on the Bethe 
lattice is achieved when \i is chosen to be U/2 due to a 
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special particle hole symmetry of the impurity problem : 

d L -> d l-v d l° ~* d l-u ' f" /-a-) /it /-ff and with 

(e, -► -(e, provided fi = 17/2. 

To study these equations we use the algorithm of Ref. 
which uses zero temperature Lanczos method to com- 
pute the Green's functions of the impurity model given 
by (Q) and iterate the model until the self-consistency 
given by (Q) is achieved. Though the results to be pre- 
sented were obtained for the case where the bath was 
represented by 5 sites, several of them were checked for 
N = 7 to establish the results seen for N = 5. The typi- 
cal number of iterations varied between 50 — 65 for robust 
convergence. We found substantially different behavior 
in the cases of strong and weak correlations (induced by a 
quenching of antiferromagnetism mimicked by ferromag- 
netic interactions in this paper) which we discuss sepa- 
rately. 



A. Strong Magnetic Correlations 

For small U, the system is in the paramagnetic phase. 
As U increases beyond a critical U cm , antiferromagnetic 
moments develop and the up and down spectral func- 
tions /Oo-(w) = (— l/7r)IraG (r (cj) are no longer equal. The 
spin up spectral function has more spectral weight in the 
upper Hubbard band than in the lower Hubbard band 
and vice versa for the down spin spectral function. The 
low frequency peak in p a (u}) is no longer centered around 
w = but is split into two peaks centered around some 
±wo with a minimum at u) = 0. This can be attributed to 
the fact that the effective staggered magnetic field that 
is generated when antiferromagnetism develops splits the 
quasiparticle bands. As U is increased further, the den- 
sity of states at u> = decreases until a critical Umit 
where p a (y> = 0) = and the system becomes insulat- 
ing. For U > Umit, a gap A which grows with U, opens 
continuously. 

A plot of the staggered magnetization m versus U is 
shown in Fig. ^. U is in units of t in all the graphs, m 
increases monotonically with U and does not exhibit any 
special feature at the MIT. Note that though the mag- 
netic moment is quite large at the MIT and increases with 
increasing tijt\ it saturates only when one is well into the 
insulating phase. In the presence of magnetic long range 
order, the self energy does not show any anomalous be- 
havior across the MIT. This should be contrasted with 
the paramagnetic case where ImS diverges as (iu> n ) 
at the MIT. To obtain an insight into the nature of the 
MIT, we use the following low energy parametrization of 
Si and £2 m conjunction with the numerical results 

Ei(kj„) = h+ (1 - -)iu n (6) 

z 

E2(kj„) = -h + (1 - -)iu n 



where h is the staggered field generated by the interac- 
tions in the antiferromagnetic phase, z goes to zero at 
the PM-PI transition, but in the presence of magnetism 
it is non-monotonic and remains non-zero even in the 
insulating phase (cf. Figj^). 




FIG. 1. The magnetization m (bold line) and z (dashed 
line) versus U for £2/^1 = 0.57 with Umit = 3.0£. 

In terms of the local Green's functions, E CT is determined 
by the following expression on the Bethe lattice 

G- 1 + J: a =tu Jn + f i a - t j-G^- t ^G a (7) 
Rewriting (0), we obtain 

1 = (tw„ + ImSi + Re£i)Gi - - |g? (8) 

1 = (iuj n + Im£ 2 + Rc£ 2 )G 2 - |giG 2 - ^G 2 2 

Using (g) in (||) and continuing to real frequencies we 
obtain the low frequency spectral functions 

Pi, 2 M =p ( 0)[l±£] (9) 

where p(0) is the density of states at the Fermi level and 
is given by 

^ = °) = dW^ 2 -A 2 (10) 

Note that p(0) does not depend on z. Since h increases 
monotonically with U, p(0) decreases with increasing U 
in the magnetically ordered phase. This yields a critical 
value h c where p(0) — * and a gap opens continuously. 
This is the point at which the MIT occurs. We, therefore, 
see that the MIT in this case is of a very different nature 
and there are no Kondo resonances which disappear dis- 
continuously at the MIT. Our results indicate that in 
the vicinity of Umit, h increases linearly with U imply- 
ing that p(0) vanishes as (Umit ~ U)? as we approach 
the transition from the metallic side. 
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We now turn to the behavior of the coefficient of the 
linear term of the specific heat 7. The basic physics con- 
trolling the behavior of the specific heat in the antifer- 
romagnetic phase is the competition between the the in- 
crease of the effective mass m* and the decrease of the 
density of states p(0). Though m* oc [1 - d ^^ ) ) = z~ l , 
initially increases, its increase is cut off by the staggered 
magnetization (cf. Fig. |l|). When the staggered magne- 
tization to is large the decrease in p(0) is the dominant 
effect, whereas the first effect dominates when the mag- 
netism is weak. The latter is true in the case of the para- 
magnetic MIT where m* and hence 7 are both oc z^ 1 
and diverge as z — > at the transition. This divergence 
is related to the fact that there is a residual entropy in 
the insulating phase which results from the spin degen- 
eracy at every site. However, such a degeneracy is lost in 
the case where the insulator is also an antiferromagnet. 
Infact using Eqs. (|6|,p|JlO|) we find 

70 z V or 

.2 

where 70 = ^f- is the specific heat coefficient of the non- 
interacting problem. Since z remains nonzero, 7 — > as 
(Umit — U)i at the MIT. As anticipated, we see in Fig.^ 
that 7 increases with U for small m and decreases for 
larger moments. 

We can generalize the above to a lattice with a realistic 
dispersion Ek = e(k) + e(k) where e and e are the contri- 
butions from the t\ and ti hoppings. The Fermi surface of 
the non-interacting system is defined by E^ F = 0. When 
the interaction U is turned on, the dispersions of the 
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FIG. 2. 7/70 vs. U for t 2 /ti = 0.57 (Umit = 3.0t). 70 
has been defined in the text. 

quasiparticles change and are given by the poles of the 
Green's function given by Eq. 2. From the preceding 
analysis, we can say that the main effect of interactions 
is to modify the band dispersion in the following manner 

E k = z[e k ±y/el + h?\ (12) 



In the non-interacting case (h = and z = 1) the two 
bands overlap and there is no gap at the Fermi surface. 
For small U in the paramagnetic phase there is no gap in 
spectrum but the bands are slightly renormalized by the 
factor z. When antifcrromagnctic moments set in, the 
non-zero h changes the band curvature and the renor- 
malized bands start moving away from each other reduc- 
ing the Fermi surface area. The manner in which the 
Fermi surface area decreases is described by the density 
of states p(0). For small h there are regions where the 
bands overlap and there is no gap in the system. At a 
critical value h = h c a gap opens up in the spectrum as 
the Fermi surface shrinks to zero signalling the MIT. 



B. Weak Magnetic Correlations: J / 



The previous scenario was characterized by relatively 
weak electron correlations, the quasiparticle residue z 
near the transition, was at most 0.65, indicating that 
a relatively high fraction of the spectral weights remain 
coherent. This can be understood in simple qualitative 
terms, when the magnetic frustration is weak, the rela- 
tively large magnetic exchange produces a large magneti- 
zation, which in turn reduces the double occupancy and 
hence the correlations. All the spin entropy is quenched 
by the spin ordering, and the presence of a hopping in 
the same sublattice favors coherent quasiparticle propa- 
gation. To access the strongly correlated regime where 
most of the spectral function is incoherent we we turn on 
a nonzero J which reduces h and m and hence increases 
Umit as suggested by ©. 

We studied the case t\ = 0.95i 2 and t\ = 0.05t 2 with 
tz/ti = 0.22. For J = 0, this system shows the usual 
band MIT discussed in the previous section, in the vicin- 
ity of U = 2t. By choosing a J(U) cf. inset of FigJ|, 
such that the magnetic moment remains very small for a 
large regime of U as shown in Fig. ||, we move the MIT 
to U = 3.8t. In the antiferromagnetic metallic phase, 
a Kondo like resonance persists in the spectral function 
up to U = 3.8i and disappears suddenly at the transi- 
tion. This resonance at small uj splits into two peaks 
when the magnetic moment becomes appreciable as was 
seen earlier. However, unlike in the paramagnetic case 
the height of this resonance is not pinned at the value of 
the non-interaction density of states but decreases with 
increasing U. We also find that the slope z decreases to 
values much smaller than was seen in the band transition 
picture as seen in Fig. |3|. z decreases until the moment 
becomes sufficiently large as to stop its reduction. In the 
insulator, which still resembles the band-insulator in that 
there is a pseudogap, the self-energy is still analytic and 
one can still define a z. 
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FIG. 3. m (bold line) and z (dashed line) vs. U for the 
case of small frustration fa/ti = 0.2 with J varying as a func- 
tion of U. Umit = 3.8t in this case. The inset shows the U 
dependence of J. 

We find that there is a jump in z and m at the transi- 
tion in addition to the discontinuous disappearance of the 
resonance mentioned earlier, reminiscent of a first order 
transition. On re-tracing the insulating solution as func- 
tion of decreasing U, we find that there is co-existence. 
The insulator survives up to U = 3.4f and there is a MIT 
of the pure band kind at that point and the solution goes 
over to one which is a metallic antifcrromagnet but with 
a much higher value of the moment m and a much larger 
z. This metallic solution disappears around U = 2At. 
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FIG. 4. 7 vs. U for the case of small frustration 
<2/ti = 0.2 with J varying as a function of U (Umit = 3.8i). 

We have measured the double occupancies in both 
these solutions and find that the metal with the smaller 
moment has higher double occupancy. This in conjunc- 
tion with the ferromagnetic coupling J ensures that the 



metal with the lower moment and higher double occu- 
pancy is the lowest energy solution. In this weak mag- 
netic correlation regime, we find that the specific heat co- 
efficient 7 is strongly enhanced cf. Fig.|]. This scenario is 
reminiscent of the behavior seen in V2O3 Note that at 
the MIT, there is a finite gap which is much smaller than 
the Hubbard gap. Also there survive coherent quasipar- 
ticle peaks albeit with very small weight. These coherent 
peaks are asymmetric about uj = in both G\ and G2. 
Close to U = A.lt these peaks vanish completely and one 
is left only with the incoherent Hubbard band structures. 
In the coexistent solution, we see that as U is reduced the 
gap closes continuously and the the weight of these two 
coherent peaks increases continuously too. 

To summarize, the character of the AM-AI transition 
is very different from the PM-PI. We found two clearly 
distinct scenarios depending on the strength of magnetic 
correlations and frustration in the system. In the limit of 
strong magnetic correlations, the transition takes place 
as a renormalizcd Slater transition. That is, up to a mul- 
tiplicative factor which remains finite at the transition, a 
gap opens continuously, and 7 — > vanish continuously 
, as the MIT is approached from the metallic side. On 
the other hand when development of antiferromagnetic 
moments in the metallic phase is suppressed a new sce- 
nario emerges: the gap from the insulating side remains 
finite at the MIT, and a substantial increase of the spe- 
cific heat is observed when the MIT is approached from 
the metallic side. Our results indicate that the transition 
is weakly first order. These two scenarios, are remark- 
ably similar to the observations reported in NiS2- x Se x 
|f|[iof and V2-J/03 || respectively and could be developed 
further to interpret other experimental observations 0j . 
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